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We present a study of inhomogeneous big bang nucleosynthesis with emphasis on transport phenom- 
ena. We combine a hydrodynamic treatment to a nuclear reaction network and compute the light 
element abundances for a range of inhomogeneity parameters. We find that shortly after annihila- 
tion of electron-positron pairs, Thomson scattering on background photons prevents the diffusion 
of the remaining electrons. Protons and multiply charged ions then tend to diffuse into opposite 
directions so that no net charge is carried. Ions with Z > 1 get enriched in the overdense regions, 
while protons diffuse out into regions of lower density. This leads to a second burst of nucleosyn- 
thesis in the overdense regions at T < 20 keV, leading to enhanched destruction of deuterium and 
lithium. We find a region in the parameter space at 2.1 x 10 -10 < 77 < 5.2 x 10~ 10 where constraints 
7 Li/H < 10~ 9 ' 7 and D/H < 10~ 4 ' 4 are satisfied simultaneously. 

PACS numbers: 26.35.+C, 98.80. Ft, 05.60.Cd 



I. INTRODUCTION 

Inhomogeneous big bang nucleosynthesis (IBBN) has 
been studied in several papers H @. In IBBN the 
baryon density is assumed to be inhomogeneous during 
nucleosynthesis. The inhomogeneity could be the result 
of a first-order phase transition occurring before BBN, or 
of some unknown physics possibly connected with infla- 
tion. 

The effects on light element production depend 
strongly on the length scale of the inhomogeneity. It 
is well known that there is a so-called "optimal scale", 
at which the production of 4 He is reduced with respect 
to the homogeneous case, due to differential diffusion of 
protons and neutrons. 

The first studies on IBBN concentrated on the reduced 
4 He production and disregarded dissipative phenomena 
other than diffusion. Later works consider also other 
transport phenomena. The collective hydrodynamic ex- 
pansion of the high-density regions was first addressed 
by Alcock et al. |19|. Jedamzik and Fuller [^(J give a 
detailed study of dissipative processes at temperatures 
ranging from T m 100 GeV to T ~ 1 keV, including dif- 
fusion, hydrodynamic expansion, and photon inflation. 

The mutual diffusion of isotopes, however, has to our 
knowledge not been properly accounted for previously. 
Diffusion of one ion species is not restricted by collisions 
with another species, if both are moving into same direc- 
tion with same fluid velocity. On the other hand, momen- 
tum transfer is enhanched between two fluid components 
flowing into opposite directions. 

Accurate treatment of transport phenomena has be- 
come increasingly important, since several estimations 
on the primordial 7 Li abundance indicate a low primor- 
dial 7 Li/H, which is difficult to accomodate in standard 
big bang nucleosynthesis (SBBN). Lithium is produced 
quite late in nucleosynthesis, and its yield is therefore 
particularly sensitive to the late-time transport phenom- 
ena such as ion diffusion and hydrodynamic expansion of 



the overdense regions. 

In this work we study inhomogeneous big bang nu- 
cleosynthesis with emphasis on ion transport. We treat 
the primordial plasma as a fluid, and handle the dissipa- 
tion of the baryon inhomogeneity through hydrodynamic 
equations. This allows us to take into account the effects 
of mutual diffusion. We discuss the hydrodynamics of the 
primordial plasma Section II. In Section III we present 
results from numerical simulations. In the last two sec- 
tions we compare the predicted isotope yields with ob- 
servations and give our conclusions. 

Throughout this paper we use the natural unit system 
where c = % = Hb = 1 • 



II. DISSIPATION OF BARYON 
INHOMOGENEITY 

A. Ions 

Consider the evolution of a density fluctuation in the 
baryonic component of the primordial plasma. We are 
interested in the temperature range T ~ 10 MeV-1 keV. 
We treat each isotope as a separate fluid. We write down 
the hydrodynamic equations for isotope i: 



drij 

dt 
dq t 
dt 



-V • (rijV 4 ) + — -| rea ,c (1) 

at 

-TVtn + F u + F ^e - n^eE (2) 



Here rii and q denote, respectively, the number and mo- 
mentum density of isotope i. We have ignored second- 
order terms in fluid velocity Vi, which is assumed to be 
small. We use the non-relativistic formula for pressure 
P = nT and assume that temperature is nearly homo- 
geneous. As pointed out in p0[ , the fluctuations in tem- 
perature are of the order of the baryon-to-photon ratio 

~io- 9 . 
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The last term in Eq. ([j]) represents production or de- 
struction of ions via nuclear rections. Terms Fi e and Fy 
represent momentum transfer due to collisions on elec- 
trons or other fluid components. The last term in Eq. 
(|J) represents an electric field, which is present if there is 
a departure from local charge neutrality. In the following 
we evaluate explicit formulae for the collision terms F. 

Scattering between non-relativistic particles. 
The momentum transfer between two non-relativistic 
fluid components close to thermal equilibrium is given by 



F kj = nrnj f f d 3 p k d 3 Pj x 



(3) 



where fk(p) is the momentum distribution of particle fc, 
such that rifc/fc(p) gives the phase space density, Uj — 
is the relative velocity, Pjk is the center-of-mass momen- 
tum, and 



'kj 



da 



(1 - cos(0))dn. 



(4) 



is the transport cross section. Assuming a small deviation 
from the Maxwellian distribution 

A-(p) = A°(p)(i + ^), (5) 

where = (uk) is the fluid velocity, we obtain 

Ffcj = -Ukn-jSkjivk - Vj). (6) 

where 



8/27>\ 1/2 . 



Here p is the reduced mass and the thermally averaged 
cross section is given by 



f 



k 2 



where k is the center-of-mass momentum. 

Let us apply the above results to neutron-proton and 
ion-ion scattering. At low energies (below a few MeV) 
the neutron-proton interaction is dominated by s-wave 
scattering. The cross section is given by p2| 



37ra? 



(a s fc) 2 + (l-|r s a s fc 2 ) 2 {a t k) 2 + (l-\r t a t k 2 ) 2 



(9) 



with a s = —23.71 fm, at = 5.432 fm, r s = 2.73 fm, and 
r t = 1.749 fm. At zero-energy limit the cross section ap- 
proaches the value 20440 mbarn. The thermally averaged 
cross section (||) can be evaluated numerically. 



The transport cross section for Coulomb scattering be- 
tween nonrelativistic ions is given by 



a* = ^{Z t Z 3 a) 2 ^A 



(10) 



where A is the Coulomb logarithm p3 24 1. The thermal 
cross section (||) becomes a 
then have 



Sn = 



Tr(Z i Z j aYA/(2T 2 ). We 
4tt f2fM ij T\ 1/2 {ZiZjCtfA 



rp2 



(11) 



Scattering on electrons. 

The collisional force excerted on a heavy particle k 
moving with velocity v/. through a thermal background 
of light particles j is given by 



f II 1 

Ffej = n k J d 3 pp J (p)^4 J (p)p = -— n fcVfc . (1 



2) 



This equation relates the force to the mobility bkj |^3|,^5| ■ 
Here pj (p) is the phase space density of particle j in the 
frame of particle fc. Assuming a thermal distribution 
p® for particle j in laboratory frame, we have Pj(p) = 
p°AE + p ■ v), and 



dp°(E) p 3 

dp — —alJp). 

V dE 3E kjKVI 



(13) 



Neutrons interact with electrons through their mag- 
netic moment. The transport cross section is H 



3tt 



a 2 K 2 



= 8.07 x 10 -4 mbarn. 



(14) 



where k = —1.91 is the anomalous magnetic moment of 
the neutron. Using MB statistics for electrons we obtain 
for the mobility 



8 / 2m e r \ 1/2 K 2 . 5 (z) 
3 V 7T J K 2 (z) ' 



(15) 



where K are modified bessel functions. 

The differential cross section for a relativistic electron 
scattering on an ion with charge Z is given by the Mott 
formula 126(1 



da (Za) 2 E 2 
dn ~ 4fc 4 sin 4 (6»/2) 



(1 -(3 2 sin 2 (6»/2)) 



(16) 



where the factor 1 — 1 sin 2 (6*/2) is the relativistic spin 
correction. This gives the transport cross-section 



a\ e (k) = MZa) 



,ml + k 2 



A. 



(17) 



where A is the Coulomb logarithm. Using again MB 
statistics we obtain 



bZ 1 = 



AttT (z 2 + 2z + 2) (Z l a) 2 Ar 
3 K2(z)e z m 2 



(18) 
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B. Electrons 

Thermal electron-positron pairs annihilate at temper- 
atures T «1 MeV-20 keV. The remaining electrons must 
be treated as one fluid component. 

For non-relativistic electrons we have 

^ = -V-Kv e ) (19) 

^ = -TVn e + F e7 + 7i e eE. (20) 

The term F e7 represents Thomson scattering on back- 
ground photons p^J , 

F e7 = -b~Jn e v e = - / ^a T t 1 n e y e , (21) 

where <tt = 665 mbarn is the Thomson cross section and 
e 7 is the energy density of background photons. 

Note that formula (|2l]) is exactly valid only well after 
electron-positron annihilation, when photon mean free 
path is large compared with the inhomogeneity scale. 
Around T 20 keV photons are still connected to the 
plasma. For precise treatment of this transition period, 
photons should be included as one fluid component. 

Ions diffusing out from the high-density regions leave 
bching a negative net charge. That gives rise to an elec- 
tric field, which forces electrons to move so as to restore 
the electrical neutrality J28| . Electrons are thus dragged 
along with ions. The motion of ions is restricted by the 
Thomson drag force (|l|), which acts on them indirectly 
through the electric field. 

If we assume spherical symmetry, the electric field at a 
given location is determined by the total charge contained 
in the spherical region closer to the symmetry center. 
The rate of change of the field is then determined by the 
flux of charge through the sphere, 




The five differential equations (§), (§), (|T|), @, and 
( p2| ) , together with the formulae for momemtum transfer 
terms, form the basis of our hydrodynamic simulations. 

C. Diffusion approximation 

It is interesting to look at how our hydrodynamic treat- 
ment relates to the common diffusion approximation, 
where the evolution of inhomogeneity is presented by a 
differential equation of the form 

— = V-(£>Vn). (23) 

Consider the steady-state solution of Eq. (|^) in ab- 
sence of an electric field, 



- TWrii - ^ niUjSij (vi - Vj) - b ie 1 n i v l = (24) 

If we ignore the motion of par ticle species other than i 
(vj « for j ^ i), equation (E4) and the continuity equa- 
tion ([l]) together lead to a diffusion equation of the form 
(p3|), with diffusion constants given by Di e = bi e T for 
scattering on electrons, and Dij = T/(rijSij) for scatter- 
ing between nuclei. 

We note here that the neutron-proton and neutron- 
electron diffusion constants calculated this way coincide 
with those given in |l8| . Also the proton-electron con- 
stant is in agreement at the limit A»l. 

The diffusion equation describes well the motion of a 
fluid if the background fluid is stationary, so that its mu- 
tual motion can be ignored. We refer to this as the ap- 
proximation of independent diffusion. This approxima- 
tion is valid in the case of diffusion of neutrons, which 
are much more mobile than the ions and electrons they 
scatter on. The diffusion equation also describes well the 
motion of ions at high temperatures (T 3> 20 keV) , where 
the dominant scattering process for ions is Coulomb scat- 
tering on background electrons. Due to their large num- 
ber the electrons can be regarded as a stationary back- 
ground. At lower temperatures the situation is more 
complicated. Ions gain or lose momentum in collisions 
on other ion components, which move with comparable 
fluid velocities. Thus the mutual motion of the fluid 
components cannot be neglected. The situation is fur- 
ther complicated by electron drag: electrons are dragged 
along with ions so that charge neutrality is maintained. 

III. SIMULATIONS 

We have written an inhomogencous nucleosynthesis 
code where a nuclear reaction network is coupled to hy- 
drodynamic equations. We assume spherical symmetry 
and use a non-uniform radial grid of 64 cells. The grid is 
adjusted according to the density profile so that the cells 
are smallest where the gradient of the baryon density 
is largest. We assume a simple initial geometry with a 
step-like density profile. The inner part of the simulation 
volume has a high baryon density rjh , and the outer part 
a low density rji. The initial conditions are determined 
by four parameters: the volume fraction /„ of the high- 
density region, density contrast R — rinj^i, the average 
density rj = f v ■ r\h + (1 — f v ) ■ rji, and the radius of the 
simulation volume r. We give r in comoving units in me- 
ters at T = 1 keV temperature. One meter at T = 1 keV 
corresponds to 4.26 x 10 6 m today. The baryon density 
is given as the baryon-to-photon ratio 77, which is related 
to fl b through r] W = 10 1O 77 = 274ft fc /i 2 . 

The code evolves 21 variables in each grid zone: the 
concentration and momentum density of e, n, p, D, 3 H, 
3 He, 4 He, 6 Li, 7 Li, and 7 Be, and the electric field. These 
variables are evolved by solving a set of 21 x 64 stiff dif- 
ferential equations in time steps proportional to the age 
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FIG. 1. Evolution of the abundances of light isotopes as 
a function of temperature, for simulation parameters r = 10 7 
m (at 1 keV), r?io = 6, f v = 0.01, and R = 10 6 . The solid 
lines show results from a complete simulation. The dashed 
lines represent a simulation where the approximation of inde- 
pendent diffusion was applied. Separation of elements leads 
to a second burst of nucleosynthesis below T < 20 keV. 



of the universe. This involves the solution of a band diag- 
onal linear system of rank 1344 at every time step. The 
nuclear reaction rates include those given in the NACRE 
compilation p9fl . The simulation is started at T = 10 
MeV and ends at T = 1 keV, or when all nuclear reac- 
tions have ceased. The final output consist of the average 
concentrations of p, D, 3 He (including 3 H), 4 He, 6 Li, and 
7 Li (including 7 Be). 

For comparison we also made a set of simulations where 
we mimicked the approximation of independent diffusion. 
We removed from the matrix all elements corresponding 
to mutual diffusion, i.e, terms that represent dependence 
of dqi I dt on . We also forced a steady-state solution. 
The electron drag was taken into account by adding to 
the momentum loss of an ion with charge Zi the Thomson 
force that would act on Zi electrons moving with the 
same velocity. This approach can be written as 



dni dni 



-TVn; 



&„•„ n,:V,; 



(25) 

b-^Z i n i w i = Q. (26) 



A. Separation of elements 

Figure |l| shows the light element abundances plotted 
against temperature, for simulation parameters r = 10 7 
m j ?7io = 6, and f v = 0.01. We compare results from 
a complete simulation, and from a simulation where we 
used the approximation (26). The complete simulation 




10' 10" 
r (m at 1 keV) 

FIG. 2. Light element abundances as a function of r, for 
simulation parameters r/m = 6, /„ = 0.01, and R — 10 6 . 
The solid lines show results from a complete simulation. The 
dashed lines represent a simulation where mutual diffusion 
was ignored. The helium mass fraction Y p (dash-dotted line) 
is shown on linear scale (right y-axis). Other isotopes are 
shown on logarithmic scale (left y-axis) as ratio of number 
density to that of hydrogen. 



shows a second burst of nucleosynthesis at T ss 20 — 10 
keV, leading to destruction of 7 Li, D, and 3 H. 

This can be understood as follows. Diffusion of elec- 
trons is inefficient at temperatures T 3> 1 keV due to 
the frequent Thomson scattering on background pho- 
tons. Thomson drag then resists the diffusion of ions 
which must drag electrons with them to maintain elec- 
trical neutrality. However, if we divide the motion of 
ions into two components, one that obeys the condition 
J2i niViZi — and thus carries no charge, and one that 
does carry charge, only the latter is resisted by Thomson 
drag. Protons and helium ions, for instance, are allowed 
to diffuse into opposite directions in such a way that the 
total charge flux vanishes. This leads to a separation 
of elements. Heavier elements tend to get enriched in 
the high-density regions, while protons diffuse out. The 
nucleosynthesis process in the high-density region is en- 
hanced by the increased concentration of heavier nuclei. 
This effect is responsible for the modified nucleosynthesis 
yields that our simulations show. 

In Fig. U we show the isotope yields as a function of 
inhomogeneity length scale, for /„ = 0.01 and rjio = 6. 
Again, we show results both for a complete simulation 
(solid lines) and for a simulation with same parame- 
ters but with the approximation of independent diffu- 
sion. The complete simulation shows a clear decrease in 
the abundances of D, 3 H, and 7 Li, as compared to the 
diffusion approximation. Also the yields of 3 He and 7 Be 
are reduced, but not as clearly. The yield of 4 He is hardly 
affected. 
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Lithium is destroyed via reaction 7 Li(p, a) 4 He. The 
mean destruction channel for 7 Be, instead, is via reac- 
tion 7 Be(n,p) 7 Li. As this reaction requires free neutrons, 
which are not available after the main phase of nucleosyn- 
thesis, 7 Be is affected little in the second nucleosynthesis 
phase. The same holds for 3 He, whose main reaction 
channel is 3 He(rt,p) 3 H. 

The reduction in D, 3 H and 7 Li due to element sepa- 
ration is most efficient at scales somewhat smaller than 
the "optimal scale" at which the 4 He yield is minimized. 
The maximal 4 He reduction occurs at a scale at which 
neutrons diffuse maximally out from the high-density re- 
gions, but the later back-diffusion is not too efficient. 
At a somewhat smaller scale, back-diffusion of neutrons 
leads to synthesis of nuclei in a narrow zone surrounding 
the high-density region. There are then plenty of nuclei 
to be transported deeper into the high-density region, 
once the separation of elements begins at T ~ 20 keV. 

Some analytic considerations may be in place here. 
Consider the steady-state solution of Eqs. (||) and (pfl|), 

- TV ii, 11,11 ,S,,:v, - Vj) - riiZieE = (27) 

3 

- TVn e -b-Jn e v e + n e eE = 0. (28) 

We are interested in the temperature regime T < 20 
keV and have ignored terms that represent scattering 
on electrons. A measure of the distance over which n e 
can deviate from y\. njZj is given by the Debye shield- 
ing distance h = (4iran e /T)- 1 / 2 @, which is orders of 
magnitude smaller than the inhomogenity length scale. 
We can therefore assume that electrical neutrality holds 
at the scale of the inhomogeneity (n e w X^ n i^)j an d 
based on that eliminate the field E. 

The evolution of the inhomogeneity is particularly sim- 
ple in two limiting cases. If the interaction between ions 
is strong compared with the electron-photon interaction 
(riiSij » b~}), as is the case at late times (T<1 keV), 
the plasma behaves as a single fluid, moving with a col- 
lective velocity v e . Taking the sum of Eqs. (E7|) and ( pgj ) 
and using the electrical neutrality we find 



T(Vn B + Vn e ) 



b e Jn e v e 



0. 



(29) 



This represents collective hydrodynamic motion of the 
plasma against Thomson drag [[l9| , [20[ . With the approx- 
imation rib ~ n e it leads to the diffusion equation for the 
baryon density, with diffusion constant Dh y d = 26 e7 T. 

In the opposite limiting case, when electron-photon 
scattering dominates over ion-ion scattering, the motion 
of electrons is suppressed by the Thomson drag. Ions 
then move under the condition that the net current car- 
ried by ions vanishes, '^2 i niZiVi = 0. Consider for sim- 
plicity a system of two ion species only, say hydrogen and 
4 He. The steady-state solution now simplifies into 

T ( Z\— — - Z 2 -— - ) = (vi - v 2 )S , i 2 (niZi + n 2 Z 2 ) 



"2 



T = 80.0 keV 



4 He i 



T = 18.0keV 



T = 6.0 keV 




T = 2.0 keV 



T = 80.0 keV 



T = 18.0 keV 



T = 6.0 keV 



X T = 2.0 keV 



(30) 
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FIG. 3. Separation of elements. Number density (left) 
and rev (right) are shown for a run with with only protons 
and 4 He present. Nuclear reaction were turned off. The sim- 
ulation started with an initial profile with R = 10, r = 10 T m, 
/„ = 0.4 3 , and uniform helium mass fraction Y p — 0.25 (up- 
permost frame). We show the density profile and flux of 4 He 
(solid line), p (dashed line), and e (dash-dotted line) along the 
radial axis of the sperical simulation volume. Helium begins 
to concentrate in the high-density region at T < 20 keV. The 
concentration of 4 He in the center reaches a maximum around 
T ~ 6 keV. The inhomogeneity is finally erased by collective 
hydrodynamic expansion against Thomson drag (last frame). 



It is now easy to see that if two isotopes have the same 
initial inhomogeneity (Vni/ni = Vn2/«2), then the one 
with smaller charge will flow into the direction of negative 
density gradient, while the one with larger charge will 
move into the opposite direction. The isotope with larger 
charge will get concentrated into the high-density region. 

In order to illustrate the separation of elements, we 
made a run with only protons and 4 He present. We 
started with a step-like initial profile with uniform he- 
lium mass fraction Y p — 0.25. The concentration profiles 
of the two elements and electrons, as well as density times 
velocity, at various temperatures are shown in Fig. |^. 
Helium begins to concentrate in the high-density region 
at T < 20 keV. At T = 6 keV the initial high-density 
region contains 62% of all helium nuclei, while in the be- 
ginning it contained 41%. The inhomogeneity is erased 
when Thomson scattering becomes inefficient in restrict- 
ing the collective motion of the plasma. 
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IV. IBBN COMPUTATIONS AND 
CONFRONTATION WITH OBSERVATIONS 

The best way to evaluate the primordial 7 Li has for 
a long time been the so-called Spite plateau H?| ob- 
served in halo stars. There is still debate on how much 
the 7 Li in Spite plateau stars has depleted from the pri- 
mordial abundance, and consequently, on the primordial 
7 Li abundance. While some authors obtain an relatively 
high upper limit 7 Li/H < 4 x 10~ 10 fl31J1 , a number au- 
thors argue for a lower value ^-Q. Ryan et al. ]35| ] 
derive the range -10.04 < log^^Li/H) < -9.72 for 
the primordial abundance. In SBBN this corresponds 
to ?7io = 10 10 ?7 < 4.2. Suzuki et al. obtain an even 
tighter range -9.97 < log 10 ( 7 Li/H) < -9.77. A recent 
study |}5| gives an intermediate estimate log 10 ( 7 Li/H) = 
-9.76 ± 0.056 ±0.06. 

The tight lithium limits of Ryan et al. and and Suzuki 
et al. arc in conflict with the low deuterium estima- 
tions from QSO absorption systems ]36|-|38|]. O'Meara 
et al. obtain the range D/H = (3.0 ± 0.4) x 10~ 5 from 
a combined study of four such systems. In SBBN this 
corresponds to 5.4 < rj W < 6.4. The tight lithium lim- 
its are also in conflict with the recent Boomerang result 
i)~6x 10- 10 |9). 

In light of the above, it is interesting to note that the 
separation of elements due to Thomson drag leads to si- 
multaneous destruction of Li and D. 

We have computed the light element abundances for a 
range of inhomogeneity parameters. Figure ^ shows the 
yields of light isotopes as a function of length scale r and 
volume fraction of the high-density region /„, for 77 = 
6 x 10~ 10 . At small scales the results converge towards 
SBBN results. The smallest 7 Li yield 7 Li/H = 10" 9 55 
was obtained at f v w lO" 15 » 0.032 and r w 7.1 m. The 
SBBN value is 7 Li/H = 10~ 9 - 39 . 

In Fig. (||) we show the isotope yields as a function of 
r and 77, for /„ = 10~ 15 and R = 10 6 . The reduction 
in 7 Li/H is more prominent at low 77, due to the fact 
that at low 77 most of the lithium is produced directly as 
7 Li, which is sensitive to the separation of elements. At 
high 77 most of the lithium comes from 7 Be, which is not 
affected as much. 

In Fig. (|^) we compare our results to a set of observa- 
tional constraints on light element abundances. The sim- 
ulation parameters are the same as in Figs. (^). For 4 He 
we choose conservative limits 0.23 < Y p < 0.25. For deu- 
terium we select constraints 10 -4 ' 8 < D/H < 10~ 44 . The 
lower limit here comes from the present D/H abundance 
in interstellar medium ^TJ. The upper limit is based 
on the two-sigma upper limit of the O'Meara estimation 
p6[ . For lithium we choose a low limit 7 Li/H < 10 -9 7 . 
We also include the constraint 3 Hc/D < 1. The lithium 
and deuterium constraints we have chosen are in conflict 
in SBBN. The upper limit to deuterium implies 7710 > 4.8 
while the upper limit to 7 Li implies 7710 < 4.2. In IBBN 
we find a region in the parameter space, where all con- 




10 10 

r (m at 1 keV) 

FIG. 6. Observational constraints on (r, 77) plane. The 
shaded region satisfies the constraints 0.23 < Y p < 0.25 (thick 
dashed lines), lO -4 ' 8 < D/H < 10~ 4 ' 4 (thick solid lines), 
7 Li/H < 10~ 9 ' 7 (thin solid line), and 3 He/D < 1 (dash-dotted 
line) . 



straints are satisfied simultaneously. The allowed region 
falls in the range 2.1 < 7710 < 5.2, corresponding to 
0.008 < Vt b h 2 < 0.019. We note that if we apply the 
approximation of independent diffusion, the allowed re- 
gion disappears. 



V. CONCLUSIONS 

We have studied inhomogeneous big bang nucleosyn- 
thesis with emphasis on transport phenomena. We com- 
bined a hydrodynamic treatment to a nucleosynthesis 
simulation. We found an effect that to our knowledge 
has been overlooked before: separation of elements due 
to Thomson drag. Thomson drag prevents the diffusion 
of the electron fluid shortly after electron-positron anni- 
hilation. Hydrogen and multiply charged elements then 
diffuse into opposite directions so that no net charge is 
carried. Helium and lithium get concentrated into high- 
density regions, which leads to enhanched nucleosynthe- 
sis and destruction of 7 Li, D, and 3 He. The effect is 
important at length scales from 10 5 to 10 9 meters at 1 
keV temperature, corresponding to 10~ 5 — 0.1 pc today. 

We computed the light element yields for a variety of 
initial inhomogeneity profiles and found a region in the 
parameter space where a low lithium constraint 7 Li/H < 



-9.7 



and a low deuterium constraint D/H < 10 are 



10 

satisfied simultaneously for 77 = (2.1 — 5.2) x 10~ 
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FIG. 4. Isotope yields as a function of the inhomogeneity length scale r and the volume fraction /„ of the high-density 
region. Other simulation parameters were r] — 6 x 10~ 1() and R — 10 6 . At the left the results converge towards SBBN yields 
D/H = 2.9 x 1(T 5 , 3 He/H = 1.0 x 10~ 5 , Y v = 0.2483, and 7 Li/H = 4.1 x 10~ 10 . 



7 




8 



Electronic address: Elina.Keihanen@helsinki.fi 

J.H. Applegate and C.J. Hogan, Phys. Rev. D 31, 3037 

(1985). 

J.H. Applegate, C.J. Hogan, and R.J. Scherrer, Phys. 
Rev. D 35, 1151 (1987). 

C. Alcock, CM. Fuller, and G.J. Mathews, Astrophys. J. 
320, 439 (1987). 

R.M. Malaney and W.A. Fowler, Astrophys. J. 333, 14 
(1988). 

N. Terasawa and K. Sato, Phys. Rev. D 39, 2893 (1989). 
H. Kurki-Suonio, R.A. Matzner, J.M. Centrella, T. Roth- 
man, and J.R. Wilson, Phys. Rev. D 38, 1091 (1988). 
H. Kurki-Suonio and R.A. Matzner, Phys. Rev. D 39, 

1046 (1989). 

H. Kurki-Suonio, R.A. Matzner, K.A. Olive, and 

D. N. Schramm, Astrophys. J. 353, 406 (1990). 

H. Kurki-Suonio and R.A. Matzner, Phys. Rev. D 42, 

1047 (1990). 

G.J. Mathews, B.S. Meyer, C.R. Alcock, and G.M. Fuller, 
Astrophys. J. 358, 36 (1990). 

G.J. Mathews, D.N. Schramm, and B.S. Meyer, Astro- 
phys. J. 404, 476 (1993). 

K. Jedamzik, G.M. Fuller, and G.J. Mathews, Astro- 
phys. J. 423, 50 (1994). 

G. J. Mathews, T. Kajino, and M. Orito, Astrophys. J. 
456, 98 (1996). 

M. Orito, T. Kajino, R.N. Boyd, and G.J. Mathews, As- 
trophys. J. 488, 515 (1997). 

H. Kurki-Suonio, K. Jedamzik, and G.J. Mathews, As- 
trophys. J. 479, 31 (1997). 

K. Kainulainen, H. Kurki-Suonio and E. Sihvola, Phys. 
Rev. D 59, 083505 (1999). 

H. Kurki-Suonio and E. Sihvola, Phys. Rev. D 63, 083508 
(2001). 

K. Jedamzik and J.B. Rehm, Phys. Rev. D 64, 023510 
(2001). 

C.R. Alcock, D.S. Dearborn, G.M. Fuller, G.J. Mathews, 
and B.S. Meyer, Phys. Rev. Lett. 64, 2607 (1990). 
K. Jedamzik and G.M. Fuller, Astrophys. J. 423, 33 
(1994). 

See, e.g., R.D. Present, Kinetic Theory of Gases, 

(McGraw-Hill Book Company, New York, 1958). 

M.A. Preston and R.K. Bhaduri, Structure of the Nucleus 

(Addison- Wesley, Reading, MA, 1975). 

L.D. Landau and E.M. Lifschitz, Physical Kinetics (Perg- 

amon, New York, 1981). 

R.J. Goldston and P.H. Rutherford, Introduction to 
Plasma Physics (IOP, Bristol, 1995). 

L.D. Landau and E.M. Lifschitz, Fluid Mechanics (Perg- 
amon, New York, 1959). 

K.R. Lang, Astrophysical formulae (Springer, 
Berlin, 1980). 

P.J.E. Peebles, Physical Cosmology (Princeton Univer- 
sity, Princeton, NJ, 1971); P.J.E. Peebles, Principles 
of Physical Cosmology (Princeton University, Princeton, 



[28] 
[29] 
[30] 

[31] 

[32] 

[33] 

[34] 

[35] 
[36] 

[37] 

[38] 

[39] 
[40] 



NJ, 1993). 

L. Spitzer, Jr., Physics of Fully Ionized Gases (Inter- 
science, New York, 1956). 

C. Angulo et al. (The NACRE collaboration), Nucl. 
Phys. A656, 3 (1999). 

M. Spite and F. Spite, Astron. Astrophys. 115, 357 
(1982); Nature 297, 483 (1982); P. Bonifacio and P. Mo- 
laro, Mon. Not. R. Astron. Soc. 285, 847 (1997). 
M.H. Pinsonneault, T.P. Walker, G. Steigman, and V.K. 
Narayanan, Astrophys. J. 527, 180 (1999); M.H. Pinson- 
neault, G. Steigm an, T.P. Walker, and V.K. Narayanan, 
astro- ph/0105439| (2001). 

S.G. Ryan, J.E. Norris, and T.C. Beers, Astrophys. J. 
523, 654 (1999). 

S.G. Ryan, T.C. Beers, K.A. Olive, B.D. Fields, and J.E. 
Norris, Astrophys. J. Lett. 530, L57 (2000). 
T.K. Suzuki, Y. Yoshii, and T.C. Beers, Astrophys. J. 
540, 99 (2000). 
P. Bonifacio et al. 



astro-ph/0204332 (2002) 



J.M. O'Meara, D. Tytler, D. Kirkman, N. Suzuki, J.X. 
Prochaska, D. Lubin, and A.M. Wolfe, Astrophys. J. 552, 
718 (2001). 

S. D'Odorico, M. Dessauges-Zavadsky, and P. Molaro, 
astro-ph/0102164 

M. Pettini and D.V. Bowen, Astrophys. J. 560, 41 
(2001). 

P. de Bernardis et al. , Astrophys. J. 564, 559 (2002). 
J.L. Linsky et al. , Astrophys. J. 402, 694 (1993); J.L. 
Linsky et al. , Astrophys. J. 451, 335 (1995). 



9 



